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Abstract 


In the present article we argue that it is possible to find numerical solution with computer 
algebra package of coupled magnetic resonance equation for describing wireless energy 
transmission, as discussed recently by Karalis (2006) and Kurs et al. (2007). The proposed 
approach may be found useful in order to understand the phenomena related to magnetic 
resonance. Further observation is of course recommended in order to refute or verify this 
proposition. 


Introduction 


In recent years there are new interests in methods to transmit energy without wire. While 
it has been known for quite a long time that this method is possible theoretically (since 
Maxwell and Hertz), until recently only a few researchers consider this method seriously. 


For instance, Karalis et al. [1] and also Kurs et al. [2] have presented such experiments and 
reported that efficiency of this method remains low. A plausible way to solve this problem 
is by better understanding of the mechanism of magnetic resonance. Meanwhile, Rangelov 
et al. [3] have presented an analogous scheme with respect to coupled-mode theory in 
quantum system. 


In the present article we argue that it is possible to find numerical solution of coupled 
magnetic resonance equation for describing wireless energy transmission, as discussed 
recently by Karalis et al.(2006) and Kurs et al. (2007). The proposed approach may be 
found useful in order to understand the phenomena related to magnetic resonance. 


Nonetheless, further observation is of course recommended in order to refute or verify this 
proposition. 


Numerical solution of coupled-magnetic resonance equation 


Recently, Kurs et al. [2] argue that it is possible to represent the physical system behind 
wireless energy transmission using coupled-mode theory, as follows: 
a,,(t)=(i@, -T,,)a,,(t)+ iKa, (0) — F,, (t) (1) 


The simplified version of equation (1) for the system of two resonance objects is given by 
Karalis et al. [1, p.2]: 
da (t) _ 





ae i(@, —il | a, (t) +ixa,(t), (2) 
mw =-i(@, — il, )a,(t)+ ika (t), (3) 

These equations can be expressed as linear 15‘ order ODE as follows: 
LO) =-iof (t) +iKg (t), (4) 
And  g'(t)=-—ißg(t)+irf (t), (5) 
Where @=(@,-iI,), (6) 
And p =(@,-—iT,). (7) 


Numerical solution of these coupled-ODE equations can be found using Maxima release 
5.23.2 as follows. First we find test when parameters (6) and (7) are set up to be 1. The 
solution is given by: 


(%01) 'diff(f(x),x)+%i*f(x)=%i*b*g(x); 
(%02) 'diff(g(x),x)+%i*g(x)=%i*b*f(x); 
(%03) desolve([%01,%02],[f(x),g(x)]); 


The solutions for f(x) and g(x) are given by: 








_ „i| [2ig (0)b + 2if (0) — 2if (0) sin(bx)] 

f(x) =e | a +f Orcos | (8) 
_ u| [2ig (0) + 2if (0)b — 2ig (0)sin(bx)] 

g(x)=e | i + Oosa) (9) 


Translated back to our equations (2) and (3), the solutions for @ = ĝ =1 are given by: 








a(t)=e" [2e (0)x+2ia, ©) = 2ia, (O)sin(xt)] 4,(0) cos J (10) 
K 
a,(t)= | Gia (0)K+2ia, Oe 2ia (OSAD a (cos, J (11) 


Now we will find numerical solution of equations (4) and (5) when a@#/#1. Using 
Maxima 5.23.2, we find: 

(%04) ‘diff(f(x),x)+%i*a*f(x)=%i*b*g(x); 

(%05)'diff(g(x),x)+%i*c*g(x)=%i*b*f(x); 

(%06) desolve([%04,%05],[f(x),g(x)]); 


And the solution is found to be quite complicated: these are formulae (12) and (13). 





[2if (O)c + 2ig(0)b — f (O)(ic +ia)] sin Je —2ac+4b° +a’) 








ic+ia)t F 
f = eto? vic? —2ac+4b? +a? (12) 





+ fO cos Ve" —2ac+4b? +a’) 





[2if (0)b + 2ig (O)a— g (0)(ic + ia)] sins Vc? —2ac+4b? +a’) 








ictia)t + 
hse ee —2ac+4b? +a’ (13) 








+ g(O)eos(5 Ve" —2ac+4b? +a’) 


Translated back these results into our equations (2) and (3), the solutions are given by (14) 
and (15), 


[2ia@,(0)c + 2ia, (0)b—a,(0) GP +ia)] sin(L é) 7 
a, (t) = e Bria) /2 z 2 +a, (0) cos &) (14) 





[2ia@, (0)c + 2ia, (0)b—a, (0) B+ ia)|sin(L é) i 
a (t) =e Pio? z +a, (0) cosl E) (15) 





where we can define a new “ratio”: 
E=Ve?—2ac+4b? +a? . (16) 





It is perhaps quite interesting to remark here that there is no “distance” effect in these 
equations, at least in theory. Nonetheless, further observation is of course recommended in 
order to refute or verify this proposition. 
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